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ABSTRACT 

Several short-lived radionuclides (SLRs) were present in the early solar sys- 
tem, some of which should have formed just prior to or soon after the solar 
system formation. Stellar nucleosynthesis has been proposed as the mechanism 
for production of SLRs in the solar system, but no appropriate stellar source has 
been found to explain the abundances of all solar system SLRs. 

In this study, we propose a faint supernova with mixing and fallback as a 
stellar source of SLRs with mean lives of <5 Myr ( 26 A1, 41 Ca, 53 Mn, and 60 Fe) 
in the solar system. In such a supernova, the inner region of the exploding star 
experiences mixing, a small fraction of mixed materials is ejected, and the rest 
undergoes fallback onto the core. The modeled SLR abundances agree well with 
their solar system abundances if mixing-fallback occurs within the C/O-burning 
layer. In some cases, the initial solar system abundances of the SLRs can be 
reproduced within a factor of 2. The dilution factor of supernova ejecta to the 
solar system materials is ~10 -4 and the time interval between the supernova 
explosion and the formation of oldest solid materials in the solar system is ~1 
Myr. If the dilution occurred due to spherically symmetric expansion, a faint 
supernova should have occurred nearby the solar system forming region in a star 
cluster. 
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1. Introduction 



The former presence of short-lived radionuclides (SLRs) in the early solar system ( 10 Be, 
26 Al, 36 C1, 41 Ca, 53 Mn, 60 Fe, 107 Pd, 129 I, and 182 Hf) has been inferred from excesses in the 
abundances of their daugh ter nuclides in meteorites, linearly correlate d with the abundance 
of a parent element (e.g., iMcKeegan fc David 120031 ; iKita et al.l 120051 ). These now-extinct 
SLRs may provide high-resolution (<0.1 million years (Myr)) chronometers for events that 
occurred during the first several million years of solar system evolution and may also be 
potential heat sources for asteroidal metamorphism and/or differentiation. 

The SLRs with relatively long mean-lives such as 107 Pd, 129 I, 182 Hf, and perhap s 53 Mn 
may have been products of steady-state nucleosynthesis in the galaxy (jJacobsen!l2005l ). while 
those with mean-lives (r) of <5 Myr, 10 Be (r=2.2 Myr), 26 Al (r=1.03 Myr), 36 C1 (r=0.43 
Myr), 41 Ca (r=0.15 Myr), 60 Fe (r=2.2 Myr), and possibly 53 Mn (r=5.3 Myr), should have 
been produced either by energetic-particle irradiation in the early solar system or by stellar 
nucleosynthesis just prior to or shortly after the birth of the solar system. Beryllium- 10, 
which was found in CAIs (calcium- aluminum-rich inclusions) (McKeegan et al., 2000), is 
not produced by stellar nucleosynthesis, but can be efficiently formed by energetic particle 
irradiation. On the other hand, 60 Fe, t he abundance of which in the early solar system 



also appears to requi re a late addition (ITachibana &: Husd 120031 ; iMostefaoui et al.l 12005 



Tachibana et al.ll2006l ). can be efficiently produced only in stars. Thus, the presence of 10 Be 



and 60 Fe in the early solar system suggests that both energetic particle irradiation and stellar 
nucleosynthesis make contributions to the inventories of solar system SLRs. However, it is 
not clear yet which process contributed more significantly to the inventory of the SLRs that 
could be synthesized by both processes, such as 26 A1, 36 C1, 41 Ca, and 53 Mn. 

There have been several attempts to find a plausible stellar source(s) for the abundances 
of SLRs in the early solar system. A low-mass (1.5 M ) thermally-pulsing asymptotic-giant- 
branch (TP-AGB) star c annot produce enough 60 Fe to matc h the initial abundance of 60 Fe 



in the solar system (e.g., iBusso et al.ll2003l ; IWasserburg et al.ll2006l ; ISahijpal fc Sonil 120061 ) . 



Models for intermediate-mass AGB stars (5 M with the solar metallicity or 3 M Q with 
1/3 x the solar metallicity) could explain the inf erred solar system abun dance of 60 Fe as 
well as the abundances of 26 Al, 41 Ca and 107 Pd ( IWasserburg et al.l 120061 ) . However, the 
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prob ability of encounters bet ween molecular clouds and AGB stars seems to be extremely 
low (IKastner fc Meyers! 1 1994 ). implying that an AGB star was an unlikely source of SLRs in 
the solar system. Mass-loss winds from Wolf-Rayet stars may have been a source of the solar 
system 26 Al, 36 C1, 41 Ca, and 107 Pd if the nuclides w ere incorporated into the solar system 
within a time-interval of ~10 5 year after production (lArnould et al.ll2006l ). However, a WR 
wind alone would not contain en ough 60 Fe and 53 Mn to explain their solar system initial 
abundances (lArnould et al.ll2006l ). 



Type II core-collapse supernovae have also been considered as plausible sources for 
SLRs. However, most supernova models imply that if a supernova provided 26 Al and 41 Ca 
into the solar system, it would also sup ply 10-100 times more 53 Mn than its estimated ini- 
tial abundance in the solar system (e.g. 



Goswami and Vanhala 2000; Ouellette et al. 2005 



Sahijpal fc Sonil 120061 ). This discrepancy could be explained by f allback of the innermost 
layers conta ining most of the 53 Mn onto a collapsing stellar core (IMeyer fc Clayton! |2000| ; 
Meyer! 120051 ) . In such a cas e, 53 Mn could primarily be derived from a different source such 
as the interstellar medium. IWasserburg et al.l (120061 ) proposed that 53 Mn and possibly 60 Fe 



were injected into the solar system as supernova ejecta with a time-interval of ~10 7 years 
after production, long enough for 41 Ca and 26 Al to decay completely, and that 26 Al and 41 Ca 
in the solar system may have been produced either by energetic particle irradiation or by an 
AGB star. 

Another problem with supernovae as sources of SLRs is overproduction of 60 Fe if all the 
26 Al in the solar system was derived from supernovae. Alth ough the yield of 60 Fe depends 
the mass loss and initial mass (e.g., iLimongi the expected amount of 60 Fe 

injected from a supernova would be, in general, a few times higher than its highest estimate 
for the early solar system, as we will show later. This problem may still remain even in 
models for fallback supernovae. 

In this study, we propose a supernova with mixing and fallback, with a kinetic en- 
ergy of explosion slightly less than that for a typical supernova (~10 51 erg) as a potential 
source of 26 Al, 41 Ca, 53 Mn, and 60 Fe in the early solar s ystem. Faint supern ovae such as 



SN1997D and SN1999br have such kinetic energies (e.g., iNomoto et al.l 120061 ) . In models 
for supernovae with mixing-fallback, the inner region of the exploding star experiences mix- 
ing, some fraction of mixed materials is ejected, and the rest undergoes fallback onto the 
core (e.g JUmeda fc Nomotolbood 120051 : INomoto et aUboodEbminaga et alil2007h . Nucle- 
osynthesis in a faint supernova with mixing-f allback successfully reproduces the elemental 



abundance patterns of hyper metal-poor stars (jUmeda fc Nomotoll2003l ; llwamoto et al.ll2005 
Nomoto et aDbood ). 
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2. Injection of supernova ejecta into the solar system materials 



The abundance of a SLR injected into preexisting solar system materials can be ex- 
press ed as follows, assuming that injected materials are well-mixed with preexisting materials 



.e.g., 



Wasserburg et al.ll2006l : ISahiipal fc Sonill2006l l: 



T 



where N R and N S1 are the numbers of SLR and a stable isotope (SI) for the initial solar 
system, respectively. Time-zero for the solar system is, in practice, considered to be the 
time of formation of CAIs, the oldest solid materials formed in the solar system. iVgjecta 
and Afgjg Cta are the numbers of SLR and SI at the time of their production, respectively. 
Nn l is the number of SI in the preexisting solar system materials. We obtained iVrHL and 



o 



cjccta 

ejecta f° r supernovae with mixing-fallback based on nucleosynthesis models by Nomoto et 
al. (2006). Their values depend on variable parameters for mixing-fallback such as the scale 
of the mixing region and an ejection efficiency of mixed materials. Details of determination of 
^ejecta an d ^f/ecta ejecta i n the context of a mixing-fallback model are shown below. N^J^^ 
and -/V^Lt c were also obtained f or non - fallback and fallback s u pernovae using nuc l eosyn thesis 



models by I Wooslev fc Weaver! (119951 ). iRauscher et al.l (12002 h 



Chieffi & Limongil (I2004h . and 



Nomoto et al.l (120061 ) for comparison. tVq 1 was taken from lAnders fc Grevessd (119891 ) . The 
variables /o an d A are free parameters in the injection model representing a dilution factor 
for supernova ejecta mixed with the preexisting materials and the time interval between the 
production of SLR and the CAI formation, respectively. 



3. Nuclides ejected from supernovae with mixing-fallback 

In the standard model for a supernova almost all the materials are ejected, but materials 
in the very inner- most region, which is typically deeper than the inco mplete Si-burning layer , 



fall back ont o the star. In the model for a supernova with fallback (IMeyer fc Clayton! 12000 



Meyer! 120051 ). such a mass-cut boundary was moved out to the incomplete Si-burning layer 
or even out to the C/O burning layer, which resulted in suppressed ejection of 53 Mn. 

In the model for a supernova with mixing-fallback, we assume two mass-cut boundaries 
at different depths of a pre-supernova star. The deeper boundary is defined as the initial 
mass-cut boundary (M cut ) corresponding to the mass-cut boundary in previous models. The 
shallower boundary is defined as the outer boundary of the mixing region (M m i x ), which is 
considered to be the volume between M cut and M mix . Most of materials in the region between 
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M cut and M m - lx fallback onto the core by gravity after complete mixing by Rayleigh- Taylor 
instabilities, but a small fraction (q) of homogeneously-mixed materials is ejected from the 
mixing region. Thus, in the present model, supernova ejecta consist of all the nuclides from 
the region outside of M mix and a fraction q of nuclides within the mixing region between 
M cnt and M mix . 

Two mass-cut boundaries (M cut and M mix ) and the fraction of materials ej ected from the 



mixing region (q) are essen tial parameters to determine N^j^ ta and iVg? ecta (e.g. Jlwamoto et al. 



2005 



cjccta. 



Nomoto et al.ll2006l ). The kinetic energies of explosion of faint supernovae (SN1997D 
and SN1999br) are estimated to be ~(0.4-0.6)xl0 51 erg, which is ~l/3 to 1/2 times smaller 
than the typical kinetic energy of supernova explosion. However, an explosi on with slightly 



Woosley fc Weaver 



less k inetic energy does not significantly affect explosive nucleosyntheses (e.g., 
19951 ). and thus we evaluated N^^. a and iV^ ecta with different sets of (M cut , M mix and q 
based on a nucleosynthesis mod el for a solar-met allic ity massive star with the kinetic energy 
of explosion of 10 51 erg (Fig. 1; iNomoto et all (120061 )). 



Regarding M cut , the initial mass-cut within the incomplete Si- burning layer suppresses 
the ejection of 53 Mn as discussed in fallback supernova models (IMeyer fc Clayton! 12000 



Meyer! 120051 ). Although a proper choice of M cut in the narrow incomplete Si-burning layer 
may explain the amounts of SLRs in the solar system including 53 Mn, very precise tuning 
of the model is required. In the present study, we varied M cut from the bottom of the 
incomplete Si-burning layer to the deeper region as in standard models to check its effect 
in the context of mixing-fallback. We found that the choice of M cut affected iVg^ a , but 
did not change significantly A^jecta as long as M cut was set deeper than the bottom of the 
incomplete Si-burning region, where no effective nucleosynthesis of SLRs discussed in this 
study takes place. As we will show below, /o is estimated to be much smaller than 1, and the 
change of iVg^ a due to change of the initial mass cut has almost no effects on the following 
discussion. We will show results for M cut placed at the boundary of the iron core, where 
peak temperatures are (10-8.5) xl0 9 K. 

The position of the outer boundary of the mixing region, M mix , was varied from within 
the incomplete Si-burning layer to the bottom of H-burning shell, and the ejection fraction, 
q, was given values between 10 _4 -10 -2 . Note that M mix and q were required to be within 
the C/O-burning region or the He-burning region and ~10 -4 , respectively, to explain the 
elemental a bundances of hyper m etal-poor stars (HE0107-5240 and HE1327-2326) by a faint 
supernova (jlwamoto et al.l 120051 ). The upper limit of q was determi ned so that the yield of 
56 Ni did not to exceed those for faint supernovae (~10~ 3 M ) (e.g.. INomoto et al.l 120061 ) . 
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4. Modeling the solar system SLRs with a supernova with mixing-fallback 

We determined fo and A to minimize deviations of modeled abundances of 26 Al, 41 Ca, 
53 Mn, and 60 Fe (see Eq.(l)) from their estimated initial abundances in the early solar system; 
26 A1/ 27 A1 = (5.0±0.5)xl(T 5 , 41 Ca/ 40 C a = (1.4±0.14)x lQ- 8 , 53 Mn/ 55 Mn = (9±4.5)xl0~ 6 , 



and 60 Fe/ 56 Fe = (7.5±2.5) xlO" 7 (e.g.. iKita et alJl2005h . Because the initial abundances of 



53 Mn and 60 Fe have not been well determined yet, we chose plausible abundances for them 
and gave larger uncertainties than for 26 Al and 41 Ca based on ranges of reported initial 
abundances. Note that the choice of plausible initial abundances for 53 Mn and 60 Fe did 
not have a significant effect on the following discussion. The uncertainties for the initial 
abundances were used as weights for minimization of fo and A. Minimization was done in 
wide varieties of fo and A, of which typical steps were 5 % and 100 years for / and A, 
respectively, and we confirmed that the grand minimum was obtained. 



5. Short-lived radionuclides from supernovae with mixing-fallback 



Calculated initial abundances of 26 Al, 41 Ca, 53 Mn, and 60 Fe normalized to meteoritic 
initial abundances for non-fallback supernovae are shown in Figure 2a-d as a function of 
the mass of the exploding star. The dilution factor, f , is estimated to be in the range of 
(0.5-6) xlO -4 , and typical time-interval, A, is ~ 0.7 Myr. As seen in previous studies, 53 Mn 
and 60 Fe are, in general, overproduced compared to other SLRs, respectively, irrespective of 
nucleosynthesis mo dels. Figure 2e shows the case for fallback models based on the nucle- 
osynthesis model of IWoosley &; Weaver! (119951 ) with no mixing, where the amounts of 53 Mn 



are much less than that expected in the early solar system (after iMeyer fc Clayton! (120001 )). 
As mentioned above, if the mass cut occurs within a narrow incomplete Si-burning layer, it 
might be possible to explain all the SLRs discussed here with fallback supernovae. However, 
it requires fine tuning of the mass cut region within the layer. Thus the general consequences 
of fallback supernova models seems to be very little ejection of 53 Mn. 

The abundances of SLRs estimated for faint supernova of different masses with mixing- 
fallback are shown in Figure 3 as a function of M mix and q. It is clearly seen that injection 
of 53 Mn and 60 Fe is suppressed and the abundances of SLRs agree with their solar system 
abundances as long as the outer boundary of mixing region is located in the C/O-burning 
layer, where peak temperatures of the shock wave range from ~4xl0 8 to ~4xl0 9 K. On 
the other hand, the modeled abundances of SLRs can not explain the s olar s ystem values if 
M mix is located in the Si-burning or He-burning layers. Ilwamoto et al.l (120051 ) proposed that 
the outer boundary of the mixing layer could be in the upper C/O-burning region or at the 
bottom of the He-burning layer to explain the elemental abundances of hyper metal-poor 
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stars, HE0107-5240 and HE1327-2326, respectively. The former case is consistent with M m - lx 
required to explain the solar system SLRs. It should be noted here again that the results 
shown here are not affected by the choice of the initial mass cut(M cut ) as long as it occurs 
in the region deeper than the incomplete Si-burning layer. 

The obtained fo and A increase from 7xl0 -5 to 2xl0 -3 and from 0.8 to 1.1 Myr, 
respectively, with increasing M mix in the C/O-burning layer. Both parameters affect the 
abundances of 26 A1 and 41 Ca, while those of 53 Mn and 60 Fe are affected mainly by /o when 
A~l Myr. Thus, f Q increases with M mix because less SLRs are ejected for larger M mix . 
However, the increase of fo to increase the amounts of 53 Mn and 60 Fe may result in the 
over-injection of 26 Al and 41 Ca, which are adjusted with increasing A in the present model. 

In addition of the calculations shown in Fig. 3, we calculated the abundances of SLRs 
with the ejection fra ction g=10~ 4 , which has been proposed for hyper-metal-poor stars 



(jlwamoto et al.l 120051 ) . but the abundance of 53 Mn was about ten times lower than the 



estimate for the initial solar system, as in the simple fallback models. 

The best estimates for the faint supernovae with mixing-fallback are shown in Figure 4, 
where the calculated abundances of SLRs reproduce their solar system abundances within a 
factor of 2. It should be emphasized here that a narrow range of mixing-fallback parameters 
(^mix and q) is not required to match the solar system abundances of SLRs. The overall 
agreement of calculated abundances with meteoritic values can be obtained for wide ranges 
of M mix and q (Fig. 3), a major advantage of this model compared to the simple fallback 
supernova models. 



6. A faint supernova with mixing-fallback as a source of solar system SLRs 

Previous supernova models, which attempted to explain the initial abundances of solar 
system SLRs, had the problems of over-injection of 53 Mn and 60 Fe compared to 26 Al and 
41 Ca. Such problems can be solved in the faint supernova model with mixing-fallback with 
the following reasons. 

Manganese-53 is synthesized mainly in the incomplete Si-burning layer, while 26 A1, 41 Ca, 
and 60 Fe are produced more abundantly in outer regions (Fig. 1). Although a small fraction 
of the 53 Mn is ejected, most falls back only to the collapsing core when the outer region of 
the mixing layer is located in the C/O-burning layer. On the other hand, some fractions 
of other SLRs, produced outside of the mixing layer, are ejected without experiencing the 
mixing-fallback. Together, these features suppress the amount of 53 Mn injected into the solar 
system materials compared to other SLRs. 
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Aluminum-26 is the only nuclide, among SLRs considered in this study, that forms 
abundantly in the He and H layers and the explosive O-burning layer of a massive star with 
the solar metallicity. If the oute r boundary of the mixing region is located in the C/O layer, 
more than 80 % of 26 Al should have its origin in the He and H layers. This is be cause 26 Al 
in the C/O layer was destroyed by the He burning (e.g., iLimongi fc Chieffil 120061 ) and most 
of 26 Al in the explosive O-buring layer falls back onto the central remnant. The amount of 
26 Al in the supernova ejecta is thus less affected by the mixing-fallback process, which lowers 
the 60 Fe/ 26 Al ratio. Aluminum-26 may be ejected in a mass-loss wind before the supernova 
explosion, which ma y explain the possible time interval between 60 Fe- and 26 Al-injection 
faizzarro et al.ll2007h . 



We conclude that 26 Al, 41 Ca, 53 Mn, and 60 Fe in the early solar system could have been 
brought from a single massive star that experienced a less-energetic supernova explosion. 
The estimated dilution factor, /q, may pro yide a geometric c onstr aint on the supernova and 
the solar system materials as discussed by ISahijpal fc Sonil (120061 ). Assuming the spherical 
symmetric ejection from the supernova, fo can be basically expressed by an injection effi- 
ciency of ejecta into the solar system materials (a) and a solid angle (O) of the solar system 
materials (a proto-solar molecular cloud core or a proto solar system disk) on a sphere of 
ejecta with a radius corresponding to the distance (D) between the supernova and the solar 
system materials. If ej ecta are injected into t he molecular cloud core (0.1 pc in radius, 1 M 
in mass, and a of 0.1 (IVanhara &: Bossll2002l )). D is esti mated to be ~2-5 pc. I n the case of 
injection into the proto solar system disk as discussed in lOuellette et al. J2005h (100 AU in 
radius, 0.01 M in mass, and a of 1), D of ~0.3-0.8 pc was obtained. 

It is difficult to evaluate which case is more plausible for the solar system forming envi- 
ronment at present. However, in either case, the supernova explosion should have occurred 
nearby the solar system materials, which supports the idea of the birth of the solar system 
in a star- cluster containing massive stars. The lifetime of the star cluster is several to 10 



Myr (e.g.. lAdams fc Laughlinll200ll ). within which only ma ssive stars (> 20-25 M ) explode 



The proportion of massive stars in a cluster is low (e.g., iKroupal l200ll ). and thus it may 
be highly unlikely that multiple supernovae brought SLRs into the solar system materials 
within a lifetime of the star cluster. In the model for a faint supernova with mixing-fallback, 
solar system SLRs with mean lives of <5 Myr can be from a single supernova. 
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Fig. 1. — Mass fractions of SLRs as a function of Lagrangian mass coordinate after the 
explosion of a solar- metallicity 40 M star with the kinetic energy of explosion of 10 51 erg 
jNomoto et al.lhood ). 
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Fig. 2. — Calculated initial abundances of 26 Al, 41 Ca, 53 Mn, and 60 Fe for non-fallback su- 
pernovae (a-d) and fallback supernovae (e), normalized to their estimated initial abun- 
dances in the solar system, as a function of the stellar m ass, (a) Nucleosynthesis m odels 
for non-fallback superno vae b y IWooslev fc Weaver! (119951 ). (6) iRauscher et al.l (120021 ). (c) 
ChiefE &: Limongil (120041 ) . (d) iNomoto et a l. (2006). (e) N ucleosynthesis models for fallback 
supernovae (based on results of IWooslev fc Weaver! (119951 )). Details of nucleosynthesis mod- 
els such as the initial mass cut and the kinetic energy of explosion can be found in references 
listed above. Error bars represent uncertainties of initial abundances of SLRs estimated from 
meteorites (see text). 
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Fig. 3. — Calculated initial abundances of 26 Al, 41 Ca, 53 Mn, and 60 Fe for supernovae with 
mixing-fallback having various masses as a function of M mix and q. The calculated abun- 
dances are normalized to their estimated initial abundances in the solar system. The dilution 
factor Jo ranges from 7xlCT 5 to 2xl0~ 3 , and time interval A ranges from 0.8 to 1.1 Myr. 
Each panel displays the results of 18-20 sets of calculations and the horizontal axis gives 
the mass of material in solar masses inside the outer boundary of the mixing region for 
each calculation. Dotted and dashed lines show the boundaries between the Si-burning and 
C/O-burning layers and between the C/O-burning and the He-burning layers respectively. 
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Fig. 4. — The best estimates for the abundances of 26 Al, 41 Ca, 53 Mn, and 60 Fe in models 
for faint supernovae with mixing-fallback. The values for the parameters are: g=0.001, 
M mix =3.1 M (6.7xl0 8 K; peak temperature of the shock wave), /„=1.90x 10 -3 , and A=1.07 
Myr for 20 M , g=0.005, M mix =3.0 M (1.4xl0 9 K), / =l-34xl0- 4 , and A=0.83 Myr for 
25 M , ?=0.001, M mix =6.7 M (7.8xl0 8 K), /o=4.35xlO" 4 , and A=0.87 Myr for 30 M , 
and g=0.005, M mix =11.8 M (5.8xl0 8 K), / =1.64xl0" 4 , and A=0.75 Myr for 40 M . 
Error bars represent uncertainties of initial abundances of SLRs estimated from meteorites 
(see text). 



